Shear banding, aging and noise dynamics in soft glassy materials 
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The 'soft glassy rheology' (SGR) model gives an appealing account of the flow of nonergodic soft 
materials in terms of the local yield dynamics of mesoscopic elements. Newtonian, power-law, and 
yield-stress fluid regimes arise on varying a 'noise temperature', x. Here we extend the model, to 
capture the idea that the noise is largely caused by yield itself. The extended model can account 
for the viscosity-bifurcation and shear-banding effects reported recently in a wide range of soft 
materials. A variant model may shed light on shear banding and strain-rate hysteresis seen in 
glassy star polymers solutions. 

PACS numbers: 



Nonergodic materials, such as colloidal gels and dense 
emulsions, have unusual rheology that is exploited in 
many commercial applications such as paints, foodstuffs, 
and cleaning products [l|. Understanding their flow 
behavior represents a profound challenge to statistical 
physics: this is intractable in its full generality, but 
promising progress has recently been made at three lev- 
els. Fully microscopic work - based for instance on 
mode-coupling theory (MCT) - is so far limited mainly 
to monodisperse colloidal glasses Q; it remains techni- 
cally formidable. Purely phenomenological approaches, 
in which a continuum stress or strain representation of 
the local state is supplemented by just one or two vari- 
ables (e.g. 'fluidity') [ail are much simpler but of rather 
limited predictive power, although some recent develop- 
ments create a more precise link between microscopic 
physics and the structural variables that control the con- 
tinuum behaviour. (An example is the recent use of a 
timc-cvolving non-affinc parameter in a model of star- 
polymer rheology [j|.) The middle ground is represented 
by mesoscopic models [l|, 0, @|, in which explicit but 
coarse-grained degrees of freedom obey some specified 
dynamics. 

Among mesoscopic approaches, the SGR model Q 
has become widely used to interpret experiments on 
soft matter, biophysics, and even conventional glasses 
[3, Hi . The qualitative successes of SGR include the pre- 
diction of power-law (er ~ -y^ -1 ) and Hcrschcl-Bulklcy 
(a — ay ~ j 1 ~ x ) regimes for steady-state flow curves, 
where a is shear stress, ay a dynamic yield stress and 7 
strain rate. These regimes are controlled by an effective 
noise temperature, x, which governs the jump rate out 
of local traps. In the yield stress regime (x < 1) SGR 
also predicts rheological aging, primarily via a so called 
'simple aging' scenario [i ljioj ]. This is something that 
MCT cannot yet capture [2j. On the other hand, an im- 
portant drawback of SGR is that its mesoscopic physics 
(trap-hopping) is somewhat generic, with few clues as 



to how the parameters of the model should be varied to 
address different classes of physical systems - for exam- 
ple, hard-sphere versus soft-sphere interactions. (Purely 
phenomenological models also suffer from this drawback, 
whereas MCT takes as input the equilibrium structure 
factor which can account for such differences.) Another 
drawback of SGR is that it does not readily admit non- 
monotonic flow curves of the type that would lead to 
shear banding - the coexistence at fixed a of layers with 
different 7, or vice versa [HI ]. Such limitations may be 
linked to the fact that the SGR model takes a; as a con- 
stant model parameter, while in practice the noise tem- 
perature for a given trap should depend on the level of 
jump activity in its vicinity. 

Currently we are not in a position to derive the x dy- 
namics from first principles, nor even mesoscopic ones. 
Rather, in what follows we treat the x evolution at a 
phenomenological level (rather as fluidity is treated in 
[H, H|)- We will show that even a very simple choice can 
then account for the 'viscosity bifurcation' seen in many 
soft materials [1, [13]. On applying a step stress, there is 
a stress threshold below which the system remains solid, 
but above which it flows homogeneously at a shear rate 
that exceeds a finite minimum value; we find moreover 
that the critical stress for the viscosity bifurcation itself 
depends on sample history prior to the application of step 
stress. At imposed shear rate, the same physics results 
in coexistence of a rigid, aging, 'cold' band (x = x± < 1) 
and a fluid, ergodic, 'hot' band (x — xi > 1) at a common 
stress a obeying = <Jy(x2) < a < cjy(xi) [l3|. Within 
a different variant of our model, we also find below a 
distinctive hysteretic shear-banding scenario reminiscent 
of one reported recently in star polymer solutions [l4| . 
Thus our work points to a possible new connection be- 
tween glass-based descriptions of nonergodic matter and 
two major experimental scenarios, complementingprevi- 
ous, purely phenomenological, modelling for these [3, EH- 
However, our approach remains semi-phenomenological 
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in that the choice of model variant and/or parameters 
cannot be linked directly to microstructure. 

The SGR model starts from a trap dynamics for meso- 
scopic elements. The jump rate out of a given trap is 
To exp[— (E — k£ 2 /2)/x], with Tq an attempt rate, E the 
well depth, k a (uniform) elastic constant, and t the lo- 
cal strain. The latter evolves between jumps as I = 7 
representing afhne shearing by macroscopic strain. New 
traps have (for simplicity) £ = 0; they have E values 
drawn from a prior distribution p(E) oc e~ E whose form 
is chosen, following Bouchaud [1 01 ] , to engineer an arrest 
transition, which occurs at x = 1. As stated previously x 
is viewed as an effective noise temperature [7|. However 
a clear physical interpretation of this parameter remains 
lacking. One idea is that a thermal nonergodic system 
will, on cooling, hover in a 'marginal state' with an effec- 
tive temperature close to that of its glass transition [lj| ; 
another is that x is mechanical noise created by plastic 
rearrangements elsewhere in the system. The latter im- 
plies a coupling between x and the flow dynamics, which 
was acknowledged but neglected in previous work (except 
fill]). an d whose consequences we now explore. 

To do so, we start with the SGR equation for the trap 
probability distribution P(E, £, t) 

P = -7^ - T e-( E - ke2 /^P + Tp{E)5(l) (1) 

(with T = T Q (e- ( - E - M2 / 2 ')/ x )p the total jump rate) and 
couple this to a relaxation-diffusion dynamics for x: 

d 2 x 

T x x(y) = -x{y) + x +S + \ 2 -^=0 (2) 

We have assumed for simplicity that x and 7 depend on a 
single spatial coordinate y in the shear gradient direction. 
The source term S{y), which represents the pumping of 
noise by jump events, clearly depends on P(E 7 £,t) at 
position y, which we henceforth denote P(y)- 

In what follows we explore two model variants, distin- 
guished according to different choices for the dependence 
of S upon P(y). Model 1 has 

S(y) = a(£ 2 /r) P{y} (3) 

where r = exp[(T5 — k£ 2 /2)/x] is the (dimensionless) trap 
lifetime. In this model the noise is pumped by dissipation 
of elastic energy. Model 2 instead has S(y) = a(l/"r)p( y ) 
such that all jumps contribute equally to noise regard- 
less of the local strain released. As discussed below, the 
underlying constitutive curve of both models has a ver- 
tical branch at 7 = 0, terminating in a yield stress; and 
a fluid branch that persists to stresses below this yield 
value (Figs. [T] and HJ. Both models are thereby capable 
of capturing a coexistence of unsheared and fluid bands 
at a common shear stress. Indeed, many of the flow phe- 
nomena that are for concreteness discussed below in the 
context of model 1 can in principle also arise in model 



2. Beyond these common phenomena, model 2 addition- 
ally allows the fluid branch to persist right down to the 
origin, allowing it to capture, for instance, the large hys- 
teresis loops seen during shear rate sweeps in glassy star 
polymers (Fig. [1] below). 

In the second equality of we have for simplicity set 
t x —> 0, so that local jump rates adapt rapidly to changes 
in nearby activity levels. In the same spirit, we neglect 
fluid inertia and therefore impose force balance, which, 
for planar Couette flow, requires uniformity of the shear 
stress: a(y) = k(£)p( y ) = a, a constant. 

The diffusive term in ([2]) represents the nonlocal ef- 
fect of jumps. More generally one might replace S in 
© by / dy'dt'S{y' ' ,t')Q{y -y',t - f) where Q is a ker- 
nel. We assume that Q is effectively of short range: most 
of the noise comes from nearby jumps. (Note that, de- 
spite the long-range nature of elastic interactions, the 
stress field caused by a randomly signed sum of plas- 
tic strains distributed through space is likewise domi- 
nated by local contributions If one further assumes 
that a smoothly varying mean activity level governs the 
rates of individual stochastic events, a gradient expan- 
sion in activity is thereby justified. To obtain ([2]) as 
written, we take the leading order correction within this 
expansion, S — > S + X 2 d 2 S/dy 2 , and iterate once to set 
d 2 S/dy 2 — > d 2 x/dy 2 . The latter replacement should be 
harmless, and greatly simplifies the numerical analysis. 

To solve the above systems numerically we discretize 
into i = l...n streamlines equally spaced in the y direction 
with a spacing A. We take on each streamline a separate 
ensemble of j = l...m SGR elements with barrier heights 
Eij. The stress on streamline i is cr, = (fc / m) J, • la . 
Periodic boundary conditions are used, and force balance 
is imposed as follows. Suppose a jump occurs at element 
ij when its local strain is I — I. By updating all elements 
on the same streamline as £ — > t + l/m, force balance 
is maintained, but with a stress level that is unchanged 
by the jump. Further updating all elements throughout 
the system as I — > £ — I J ran restores the global stress 
to the proper reduced level. This algorithm can also be 
thought of as the r\ s — > limit of one in which a small 
Newtonian viscosity r\ s is introduced alongside the elastic 
stress of SGR elements, and local strain rates are set to 
maintain a(y) + 77^7 = a at all times. The SGR sector of 
the numerics is handled by a waiting-time Monte Carlo 
(WTMC) algorithm [HI that chooses stochastically both 
the element and timing of the next yield event. After 
each such event, force balance is applied; every N WTMC 
steps, is then evolved to steady state (we found TV ~ 
10 sufficient for accuracy, and more efficient than TV = 1). 

We now discuss our results for Model 1, letting To = 
k = 1 define time, strain and stress units. First, note 
that in a homogeneous steady state ([2|) is equivalent to 

x — x + 2aa(x, 7)7 (4) 

which allows us to generate constitutive curves cr( A f) 
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FIG. 1: Dashed line: constitutive curve 17(7) for xo = 0.3, a = 
2.0. Symbols: WTMC results for quasi-stationary stress a at 
various imposed mean strain rates 7 found after start up of 
steady shear. Shear banding is present whenever the WTMC 
data deviates significantly from the dashed curve. Dotted 
lines: viscosity bifurcation points (±0.05) in step stress for 
log 10 t w = 1.0, 1.5, 2.0, 2.5 (bottom to top); c/FigEJ 



from the SGR constitutive equation [7( without recourse 
to WTMC. The result is shown in FigQ] for the case 
xq = 0.3 and a = 2.0. Also shown is our WTMC data 
for the actual flow behavior in planar Couette flow. This 
is quasi-stationary data, taken during a strain window 
150.0 < 7 < 200.0 following startup of steady shear. 
As in any strain-controlled experiment, only the mean 
shear rate 7 is imposed: the system is free to choose a 
banded state, and indeed it does so for a wide window 
of strain rates (Figs ll|2| ). All these data were generated 
with n = 100, to = 1000 and A/A = 0.5 or 1.0 (suffi- 
cient to resolve the interface). For each run the sample 
is 'fresh', quenched from a state with P(E,£) ~ p(E)S(£) 
at t = t w ~ 1; waiting for a time t w 3> 1 at x = x$ before 
startup of shear yields a larger stress overshoot but the 
same plateau stress. Note that the multiple bands seen 
at 7 = 0.2 in Fig. [2]can only arise in a planar shear flow: 
they would presumably be eliminated by any small cur- 
vature in the flow geometry, as is almost always present 
exp er imcnt ally. 

The results in FigQ] are consistent with the following 
picture. The homogeneous flow curve has a yield stress 
a y = uy(xq) inherited from the simple SGR model. 
Unlike simple SGR, however, 17(7 > 0) is thereafter a 
decreasing function, before re-stabilizing at higher shear 
rates. (In fact the initial slope of the curve is positive, 
but for the chosen parameters (7(7) has a maximum at 
a — ay ~ 2 x 10~ 4 , invisible in Fig[TJ and falls back below 
<jy for 7 > 10~ 3 .) This creates the standard conditions 
for shear banding but, unusually, the more viscous band 
has 7 = This band is thus effectively solid, with 
local strain rate that is close to zero (Figj2|), and slowly 
decreasing with time. The latter represents an aging ef- 
fect: this band has a low noise temperature x ~ xq = 0.3 
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FIG. 2: Profiles of noise temperature x(y) and velocity v(y) 
for xo = 0.3 and a = 2.0 with imposed mean shear rate 7 = 
0.05,0.1,0.2 (dotted, dashed, solid). Data is time averaged 
over the second half of a run with total strain 100. 

such that <Jy(x) > c, which is the condition for noner- 
godicity to arise within SGR. This aging can be confirmed 
directly by study of the correlator C(t,t w ) (the fraction 
of unhopped particles) within the low shear band. This 
shows simple aging behavior, as does the standard SGR 
model at the given x = xq In contrast, the high-shear 
band has high activity, self-consistently maintaining it in 
an ergodic state of high x and low viscosity. 

This scenario is supported by simulations on molecu- 
lar models fl3| . and is consistent with experiments on 
viscosity-bifurcating materials [!, EH. Indeed, we show 
in Fig|3] a plot of viscosity n = <r/j against time t fol- 
lowing a step stress on fresh (t w = 1) samples. As 
expected, without fine-tuning of a, we find flows that 
are homogeneous. There is a clear bifurcation between 
a < <t c , for which n(t) is large and increases without 
limit with t, and a > a c for which r](t) falls onto a 
plateau. Close to a c the time to nucleate the non-flowing 
branch is long for large systems; but assuming this re- 
mains finite our best numerical estimate for fresh samples 
is <j c j — 0.675 ± 0.005. Intriguingly, this is very close to 
the shear-banding stress plateau found in FigJT] and well 
below the dynamic yield stress o~y(xq) for homogeneous 
flow. For fresh samples, the strain-rate overshoot follow- 
ing step stress at a c j < a < cry carries the system onto 
the fluid branch, where it remains. 

The inset to Figj3] shows another intriguing facet of 
the viscosity bifurcation: the critical stress is itself an in- 
creasing function of t w (see also FigJT]). Thus, if subjected 
to step stress a c j < a < uy{xo), a scries of samples of 
different ages will show a bifurcation between flow (small 
t w ) and arrest (large t w ). In an SGR context this is nat- 
ural: in an aged sample, elements have fallen into deeper 
traps, so that the strain-rate overshoot is more modest, 
and the transient yield events caused by the step stress do 
not receive enough feedback to cause runaway to a fluid 
state. Although not emphasized in 0], a similar t w bi- 
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FIG. 3: Viscosity bifurcation for xo = 0.3, a — 
2.0, n = 100, m = 1000. Main figure: age before shear 
t w — 1; stress values (top to bottom at right) a — 
0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.64, 0.65, 0.66, 0, 67, 0.68, 0.69, 0.70, 
0.71,0.72,0.73,0.74,0.75,0.76,0.77,0.8. Inset: stress a = 
0.72, age before shear log 10 t«, = 0,0.5,1,1.5,2,2.5,3,3.5,4. 
Dotted line shows the asymptote calculated semi-analytically. 




FIG. 4: Dashed lines: constitutive curve (7(7) for Model 2 
with x = 0.15, a = 3.75. Symbols: WTMC data (n = 
100, m = 1000, A = 0.5A) for up/down strain-rate sweeps 
(left-pointing and right-pointing triangles respectively); t r = 
200,400,800 (thin, medium, bold symbols), initialized in a 
homogeneously aged state of t w = 10 4 . Dotted line, as guide 
to the eye, the quasi-steady stress attained at long times in 
shear startup for 7 < 0.1 in the shear banding regime. 



furcation holds for the phcnomcnological model reported 
there, in which an inverse fluidity or 'jamming param- 
eter', A, evolves with time. Our SGR results broadly 
support the idea Q that this evolution represents aging. 

We next turn to Model 2, focussing on a specific regime 
that may be relevant to experiments on star polymers 
which showed, under relatively rapid strain-rate sweep 
(residence time t r = 10s per observation point), an ap- 
parently conventional, monotonic flow curve. However, 
for slower sweeps (t r = 10 4 s) a much larger, almost con- 
stant, stress was measured at small 7. Shear banding 
in this region was confirmed by NMR vclocimetry, with 
7 ~ in the slow band [l5| . A strong hysteresis was also 
seen, the less viscous branch persis ting to much lower 
strain rates on the downward sweep [14| . 

Fig.Hshows data for Model 2 with x Q = 0.15, a = 3.75. 
The right-pointing triangles show the stress response to 
a slow upward strain rate sweep for a sample that was 
in a homogeneously aged state of t w = 10 4 before shear. 
A stress plateau is clearly seen for shear rates 7 < 0.3. 
In this regime, the system forms coexisting glassy and 
Sowing shear bands. At the lowest applied shear rates, 
the stress does not have time fully to attain the plateau 
value before the strain rate is swept on to a higher value. 
This accounts for the reduced stresses at the far left hand 
edge of the plot. As would be anticipated, this reduction 
is pushed to smaller shear rates for slower ramp rates. 
For shear rates 7 > 0.3 the system flows homogeneously 
on the fluid branch of the constitutive curve. 

A remarkable feature of Model 2, not seen in Model 1, 
is that the constitutive curve <j(j) remains multivalued 
down to 7 = where a quiescent glass at x = xq and a 
fluid phase at x% = xq + aT(xi)/To, both exist [19(. If a 



system is prepared on the fluid branch at high 7, then 
barring macroscopic nucleation events (which do not oc- 
cur at measurable rate in our numerics), the shear rate 
can be ramped down to zero while maintaining homoge- 
neous fluidity. The same applies to the downward part 
of an upward then downward shear rate ramp, as shown 
by the left-pointing triangles in Fig. 0] Apart from the 
absence of nucleation, which we assume would cause a 
swecp-rate-dependent escape from the fluid branch on 
down-ramping the strain-rate, the WTMC data in FigH 
intriguingly resembles that reported experimentally in 

Note more generally that nucleation events are hard 
to capture in numerics (see e.g. [20]) and their absence 
from Model 1 might give, for instance, small shifts in the 
independent a c found there. However, the qualitative 
features of both models appear robust to such effects. 
Numerical difficulties also impede detailed study of the 
limit Xo — > 0, which might be the realistic limit when 
true thermal noise is negligible 0- While qualitatively 
similar flow curves should arise, this limit of the standard 
SGR model is semi-deterministic (the element of lowest 
barrier height is always next to yield - reminiscent of 
extremal models of self-organized criticality) so that the 
treatment of yield-induced noise becomes more delicate. 

In conclusion, we have presented an extended SGR 
model for soft glasses under flow, in which the noise tem- 
perature x varies to reflect the dependence of local jump 
rates on yield events elsewhere within the sytem. Two 
versions of the model offer a connection to the physics 
of viscosity bifurcations (Model 1) and star polymer col- 
loids (Model 2), neither of which could previously be ac- 
counted for within a simple (uniform x) SGR framework. 
So far, our numerical work has addressed only the case of 
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periodic boundary conditions. Very recent experiments 
however emphasise the perturbing influence of bound- 



aries 



21], which might or might not be linked to their 



effects on the noise temperature; we hope to explore this 
in future work. Meanwhile, we believe that mesoscopic 
models, such as those developed here, continue to offer a 
useful compromise between first pri ncip les and fully 
phenomenological approaches [3, 13, 14 1 to the rheology 
of nonergodic soft matter. 
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